EQUATIONS OF CAMASSA-HOLM TYPE AND JACOBI 
ELLIPSOIDAL COORDINATES 



K.L. VANINSKY 

Abstract. We consider the integrable Camassa-Holm equation on the line 
with positive initial data rapidly decaying at infinity. On such phase space we 
construct a one parameter family of integrable hierarchies which preserves the 
mixed spectrum of the associated string spectral problem. This family includes 
the CH hierarchy. We demonstrate that the constructed flows can be interpreted 
as Hamiltonian flows on the space of Weyl functions of the associated string spec- 
tral problem. The corresponding Poisson bracket is the Atiyah-Hitchin bracket. 
Using an infinite dimensional version of the Jacobi ellipsoidal coordinates we 
obtain a one parameter family of canonical coordinates linearizing the flows. 



1. Introduction. 

1.1. Arnold's problem. Separation of variables is the simplest and the most 
powerful integration method for equations of motion in classical mechanics. Of 
course, there is no a general rule which allows one to find such separating coordi- 
nates in the general setting. "Therefore we have to go in the opposite direction and 
knowing some remarkable substitution to find a problem where it can be success- 
fully applied", [12]. One such remarkable substitution, known as Jacobi ellipsoidal 
coordinates is widespread. Below is a list (incomplete) of classical problems where 
this substitution can be applied, [2]. 

• Plane motion in the field of two attracting centers (Euler 1760). 

• Kepler's problem in the homogeneous force field (Lagrange 1766). 

• Geodesic motion on n-dimensional ellipsoid (Jacobi 1866). 

• Motion of uncoupled harmonic oscillators constrained to move on the n- 
dimensional sphere by the field of a quadratic potential (Neumann 1859). 

The following question 1 was posed by V.I. Arnold: 
Problem. 1981-29. [1]. Generalize Jacobi ellipsoidal coordinates to the infinite 
dimensional setting. Find equations of mathematical physics integrable by this 
method. 

The goal of the present paper is to demonstrate that the Camassa-Holm equation is 
an example of a PDE integrable with this technique. The Camassa-Holm equation, 
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[8], is an approximation to the Euler equation describing an ideal fluid 
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in which t > and — oo < x < oo, v = v(x,t) is velocity, and R is inverse to 
L = l-d 2 /dx 2 i.e. 



The CH equation is a Hamiltonian system m* + {to, Ti} = with Hamiltonian 



The integration of the CH equation is based on the spectral theory of an asso- 
ciated Sturm-Liouwille operator. In connection with this spectral problem Weyl, 
[28], in 1910 introduced what is now called the Weyl function. The Weyl function 
depends on the spectral parameter and defined such that some linear combination 
of fundamental solutions and its derivatives vanishes at the end of an interval. The 
novel point here is that the Weyl function carries a natural Poisson structure, the 
Atiyah-Hitchin bracket. It corresponds to the Poisson structure 1.1 on the phase 
space. The fact that Jacobi elliptic coordinates are separating coordinates for 1.1 
is proved by contour integration. 

It seems that the CH equation on the line is the simplest among all integrable 
PDE. The reason for this is a degeneration of its spectral curve. Such degenerate 
reducible rational curves appear in the compactification of the space of smooth 
hyperelliptic curves. It is important that the constructions of the present paper 
can be extended to other basic integrable models, i.e., the KdV/KP and NLS 
equations, [26, 27]. 

The rest of the introduction describes what is happening for the CH equation. 




— oo 



Introducing the function to = L[v] one writes the equation in the form 2 

to* + (mD + Dm) v = 0. 




and the bracket 




(1.1) 



2 We use notation D for the ^-derivative and • for the i-derivative. We use 5 for the Frechet 
derivative. 
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1.2. The Camassa— Holm hierarchy. We consider the CH equation with non- 
negative (m > 0) initial data and such decay at infinity that: 

" + 0O 

m(x)e^dx < oo. 

> 

We denote this class of functions by M. For such data a solution of the initial 
value problem exists for all times, see [9]. An initial profile moves from the left 
to the right and changes its shape under the flow. We also need a subclass of 
functions M.q C M. which vanish far enough to the left. Evidently A4o is invariant 
under the flow. 

The CH Hamiltonian is just one of infinitely many conserved quantities of motion 

n x 
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The first integral 
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mdx = momentum 



produces the flow of translation 

m* + {m, 7i{\ = m° + Dm = 0. 

The second is the CH Hamiltonian H2 = 7~t- The integrals TC3, TC4, etc., produce 
higher flows of the CH hierarchy. 

An interesting feature of the Camassa-Holm equation is that the invariant mani- 
fold, or in other terms the spectral class specified by these integrals, is not compact. 
As we will see this is reflected in the non-periodicity of the angle variables which 
occupy the entire real line. 

1.3. The spectral problem. Action-angle variables. One can associate to 
the CH equation an auxiliary string spectral problem, 3 

no + HOfio = 0, -2 < e < 2. 

The background information for this spectral problem can be found in [13, 15, 10]. 
The variables £ and x are related by 

^ 1 x 

x — > £ = 2 tanh — . 



3 We use prime ' to denote ^-derivative. 



Also the potential g(£) is related to m{x) by the formula </(£) = m(x) cosh 4 |. For 

initial data from M. the total mass of the associated string is finite g(Od£ < 00 • 
Evidently, for initial data from Mo there is an interval of length / where the 
potential vanishes: = 0, £ G [—2, —2 + 1]. 

The solutions of the string spectral problem are continuous, but the left f!_(£) 
and right /+(£) derivatives at some point £ may be different due to a concentrated 
mass at this point. Two solutions </?(£, A) and A) of the string spectral problem 
play an important role in the whole discussion. They are specified by the initial 
data 

V?(-2,A) =1 V(-2,A) = 

<//_(-2,A) =0 #_(-2,A) = l. 

All potentials of the string spectral problem with fixed Dirichlet spectrum 

A n : V(2,A n )=0, n = l,2,...; 

constitute a spectral class. The variables 

4 = -^-, #n = logy?(2, A n ) 

are called "action-angle" variables. The angles 9 n take values in R; they are not a 
cyclic variables. 

It turns out that the Dirichlet spectrum is preserved under the flows of the CH 
hierarchy. The flows are linearized in terms of angle variables. This follows from 
the canonical relations 

{I k , I n } =0, 
{Ok, 9n] = 0, 

{Ok, I n } = <5fc ; 
and the trace formulas of McKean, [20]: 

n n 
n n 

We note that the CH flow was linearized in terms of the so-called coupling con- 
stants by Beals, Sattinger and Smigielski, [4]. 

4 



1.4. Mixed boundary conditions. As we noted the CH flows do not change the 
Dirichlet spectrum. It is natural to consider the mixed boundary condition 



ay?(2,A) + 6^(2, A) = 0, 

where a, b are real constants. The set of all boundary conditions, i.e., the pairs 
(a, 6), constitutes MP, real projective space. For example, the Dirichlet boundary 
condition corresponds to the point (0, 1). Our goal is to associate to each bound- 
ary condition (point of MP) the family of Hamiltonian flows which preserve the 
spectrum corresponding to this boundary condition. It turns out that the space 
MP can be covered by two charts. 

The first chart corresponds to the boundary condition: 

V>(2,A)-CV(2,A) = 0, 

in which C < I. It can be constructed for all initial data from M.. The family of 
roots A n = A n (C), n — 1,2,...; produces generalized "action-angle" variables 

4(C) = "y^t; 0n(C) = logy?(2, A n (C)). 

For each value of C < I one has the family of integrals. Here are the first two 4 : 

Qt r+oo ^ r+oo 

Wi(C) = — A ~ / m(x)(l + e~ x )dx + / m(x)dx; 

4 - C J-cc 4 - C 

1 / r+oo r+oo \ 2 

m 2 (C) = ( 4 _ C y \ C J rn{x){\ + e~ x )dx - 4 J m(x)dxj + 

+00 \ 2 +00 

J dim{i) -2C J d^m(0D- 1 m(0 + 
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C-4 



+00 

+ (4 - C) f d£m(0v(0 



Evidently for C = we have classical action-angle variables and the CH Hamilto- 
nians. The Hamiltonians produce the flows 5 

m* + {m,H k (C)} = 0, k — 1,2, — 



4 The anti-derivative D 1 is denned by the formula D 1 f(x) = ^ 



+00 



-DC 



2 

5 The question of existence will be discussed at the end of the introduction. 
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These flows preserve 1(C) and move 9(C) linearly . This follows from the trace 
formulas 

n 

n 2 (c) = \Y,%{C)\ 

n 

and the canonical relations for the variables 1(C) and 9(C). 
Another chart corresponds to the boundary condition 

^(2,A) + i^>(2,A) = 0, 

in which F > —l/l. It is constructed for initial data from Mo. For F = the 
roots /Xfc(O) form the second spectrum of the string spectral problem 

fx k : ip(2,fi k ) = fc=l,2, .... 

The roots fi k (F), k = 1,2,...; produce J(F) and r(F), the second family of 
canonical coordinates: 

Jk(F) = ~^py n(F) = log^(2,/i fc (F)). 

Again, for F > —l/l one can write a family of integrals. The first two are 

1 r+°° 4p r +QO 
TAF) = / m(x)(l + e~ x )dx + / m(x)dx; 

1 / /-+oo /•+oo \ 2 

47 " 2(F) = (1 + 4F)2 [J rn(x)(l + e-*)dx + ±F J m(x)dxj - 

1 + AF 



1 + 4F 



x 



-1„ 



dxm(x) — 2 dx e x m(x)D m(x) 



+oo 

-<l+«ljf «""»<*)«(*> 

— CO 

These Hamiltonians produce the flows 

m' + {m,T fe (F)} = 0, fc = 1,2,3,.... 

The trace formulas have the form 

up) = -J2 J -(n 

n 

UF) = \j2 J n( F )> ^ 

n 



The formulas becomes especially simple when F — 0. The first integral 

f+OO 



/+oo 
m(x)(l + e~ x ) dx, 
-oo 



produces the flow 

m" + {m, 7i(0)} = m* - 2me~ :c + (1 + e~ x )Dm = 0, 

which is linearized in the variables J(0) and r(0). The flows corresponding to 
higher Hamiltonians are nonlocal. 

1.5. Ellipsoidal coordinates. The Atiyah— Hitchin bracket. The key obser- 
vation in the construction of the generalized action-angle variables / and and 
J and r is that the spectra A n (C) and Pk(F) can be interpreted as Jacobi ellip- 
soidal coordinates. This allows one to avoid any use of the hierarchy of differential 
equations as well as the string spectral problem. As a result the dynamics can be 
reformulated purely in the language of meromorphic functions on CP. First, let us 
explain the construction of I{C). 

For the string Si, i.e., for the string with free left end and fixed right end, we 
define so-called the Weyl function 

where / > is an interval free of mass, Ok > and ^2 Okj P-k < oo. The roots of the 
equation Q (X) = C, where C < I are the points A n (C). This is a classical way to 
define ellipsoidal coordinates. The construction is depicted on Figure 1. We have 
the sequence of points = Afc(C), k = 1,2, . . . which depend on the constant 
C . They interlace the poles of f2o(A): 

Hi < Ai(C) < (i 2 < A 2 (C) < . . . . 

The action variables are defined by the formula 

L(c) = - xjcy n=1 ' 2 '"" 

Classical action variables correspond to the case C = 0. 

Now we can give a similar description of the variables of J(F). We assume that 
m(x) E Mo and therefore I > 0. For the string with fixed left and right ends 
we have the Weyl function 

Eo{X) ~ "?MA) " ~WA) ~ "T + V 

where pt > and Pk/^k < oo. The points (ik{F) appear as roots of the equation 
E (X) = F, F> -I /I. This is shown on Figure 2. 
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FIGURE 1. The graph of function Q - 
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Figure 2. The graph of function E { 



We have a sequence of points Hk = f^k(F) interlacing the poles of E (\): 

m(F) < Ai < n 2 (F) < .... 
The action variables J(F) are defined by the formula 

J n {F) = 1 —, n = l,2,.... 

The second spectrum corresponds to the case F — 0. 
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The next step is to relate the functions fl and E to Poisson geometry. First 
we change the spectral parameter by the rule A — > A' = — 1/A. Note that this 
transformation maps the spectrum into the corresponding action variables. The 
Poisson bracket 1.1 for the function flo corresponding two different values of the 
spectral parameter is given by the formula 

A — jj, 

This is standard Atiyah-Hitchin bracket in the form found by Faybusovich and 
Gehtman, [11]. The bracket (or more precisely the symplectic form) was originally 
introduced for rational functions in terms of their singularities in [3]. In [11] coor- 
dinate free form of the bracket for rational functions was found. This coordinate 
free form is identically the same for an infinite-dimensional case which we consider 
here. The bracket is invariant under linear-fractional transformations. Namely, 
for any function Cl defined as 

^ _ aQ + b 
cQ + d ' 

where a, b, c and d are real constants, the Poisson bracket is given by the same 
formula, see [26]. In particular, the bracket for the function E is given by the 
same formula. The canonical relations for the generalized action-angle variables 
1(C) and 0(C) and also J(F) and r(F) are proved by contour integration. 

1.6. The spectral curve and Baker— Akhiezer function. It is instructive to 
present an algebraic-geometrical approach to the inverse problems with pure dis- 
crete spectrum developed in [18, 25]. Such problems include finite Jacobi matrices 
[18, 25], the quantum mechanical oscillator [21] and the string spectral problem. 
The algebraic-geometrical viewpoint will provide an additional geometrical insight 
into the nature of the flows constructed from the Hamiltonians 7~C(C) and T(F). 

The Riemann surface Y (Figure 3) associated with the spectral problem consists 
of two components T_ and r + , two copies of CP. A point on the curve is denoted 
by Q = (A, ±), where A e CP and the sign ± refers to the component. An infinities 
of I~± are denoted by P± correspondingly. The components are glued together at 
the points of the Dirichlet spectrum (A&, ±). 

The Baker- Akhiezer function e(£, Q) is defined on T and depends on the variable 
£, — 2 < £ < +2 as a parameter. It is defined by the formula 

un ,_\ iK&A) : if Q e T + , A = X(Q); 

el4 ' Wj " \ 4^(£,A) + ^,A), : if Qer_, A = A(Q); 

y?(2,A) 



where 

E (\) 
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Figure 3. The Riemann surface. 



The Baker- Akhiezer function has essential singularities at the infinities P± and 

simple poles at the points of the divisor 7^ = (/i^, — ), k — 1, 2, It is holomor- 

phic everywhere else. The Baker-Akhiezer function satisfies the gluing condition 

e(£,(A fc ,+)) = e(£,(A fc ,-)), k = 1,2, — 

When Q e T_ we have 

This type of formula expressing the Weyl function through the Baker-Akhiezer 
function is a starting point for consideration in the case of hyperelliptic curves, 
[26, 27]. 

The flows of the Camassa-Holm hierarchy constructed from the Hamiltonians 
Tl(0) preserve the moduli k — 1,2, ... , but move points of the divisor 7^, k = 
1,2,.. .. The flows constructed from the Hamiltonians T(0) preserve the divisor of 
the poles and move the moduli. All other flows constructed from the Hamiltonians 
H(C) and T(F) move both the moduli and the divisor. 
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1.7. The spectral class. Before going into the infinite-dimensional problem of 
the description of the spectral class of a string we present its finite-dimensional 
counterpart. Consider a finite Jacobi matrix 



v c 

C Vi Ci 








CN-3 








Vn-2 Cn-2 
C N - 2 V N -i 



c k > 0. 



All matrices with the same spectrum A < . . . < \n~i as L constitute a spectral 
class of L . There are no restrictions on the spectrum in this finite-dimensional 
case. The fact that the spectral class is diffeomorphic to R 1 ^^ 1 was already noted 
by Moser, [22]. Tomei, [24] showed that it can be compactified and becomes a 
convex polyhedron. The symplectic interpretation of this result as a version of 
the Atiyah-Guillemin-Sternberg convexity theorem was given by Bloch, Flaschka 
and Ratiu [5]. In [25] we proved that points of the divisor can serve as natural 
coordinates on the spectral class and described their range. In this paper we give 
a similar description of the string spectral class in terms of the divisor. 
First we note that the two spectra 



Ai < A 2 < A 3 < . . . ; 



(1.4) 



and 



Hi < H2 < A*3 < • • • ; (1-5) 

appear as poles of the functions Eq and Qq respectively. Their range is described 
by two theorems of Krein. The first theorem gives a description of all singular 
Riemann surfaces which might appear for initial data from Ai. 

Theorem 1.1. [16] For a given increasing sequence 1.4 to be the spectrum of a 
string So it is necessary and sufficient that there exists a finite limit 

lim ■ 

n^oo 



An 



and 



+00 ^ 



< 00, 



+00 / 

o(A )= n 

71=1 V 



1 - 



A r , 



The next Theorem gives a description of all possible divisors. 



Theorem 1.2. [16] For a given increasing sequence 1.5 to be the spectrum of a 
string Si it is necessary and sufficient that there exists a finite limit 

n 



lim 

n— >oo 
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■ w =n(i-£)- 

71=1 X ' 

Now we can present a description of the spectral class which is also due to Krein. 

Theorem 1.3. [15] For two increasing interlacing sequences of X's and /i's 

< iii < Ai < fj, 2 < ■ ■ ■ 

to be the two spectra of a string with finite mass and length it is necessary and 
sufficient that 

Now we have a complete description of the spectral curve and divisor. We 
already introduced the formal flows corresponding to the Hamiltonians 7~C{C) and 
T(F). We can say that the flows exist if the corresponding sequences 1.4 and 1.5 
satisfy the conditions of Theorem 1.3. To write an explicit form of the equations 
describing the motion of the moduli and the divisor is a routine exercise. 

Finally, we note that the problem of the description of the spectral class for the 
quantum mechanical oscillator, [21], remains open. 

1.8. The range of ellipsoidal coordinates. Now we want 

• given /i's, to describe the range of variables A(C) when C < I; 

• given A's, to describe the range of variables fi(F) when F > —l/l. 

An answer to these questions is given by two recent results. 

The first is the theorem of Kostuchenko-Stepanov, [17], a direct generalization 
of an ancient (1857) formula of Boole, [6]. In the context of the functions Eq and 
Qo their result takes the following form. 

Theorem. If YlkPk < oo in formula 1.3 or Yl/k a k < oo in formula 1.2 then the 
corresponding sequence of fik(F) or Ajt(C) satisfies 

k 

or 

^2 ~ X k(C)] < oo. 
k 

The conditions of the theorem are not only necessary but also sufficient for the 
sequence of /i's and A's to arise as coordinates for some summable sequence of p's 
or cr's. 

Unfortunately the condition of summability for the sequence of A's and /i's is 
too restrictive for our purposes. It never holds for initial data from the phase space 
M. We need the second result which due to Nazarov and Yuditskii, [23]. 
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and 



+oo 



E 



1 l/ n /2 \DUll n )\ 



< OO, 



Theorem. If Yl,Pk/^k < oo or J2<r k /fi k < oo, then 



?mV 1<0 °' 



or 



1 < oo. 



The conditions of the theorem are not only necessary but also sufficient for the 
sequence of A's and /i's to arise as coordinates for some summable sequence of 
p/A's or cr/yu's. 

1.9. Organization of the paper. In Section 2 we introduce an auxiliary Sturm- 
Liouville eigenvalue problem. Section 3 transforms this spectral problem into the 
Krein string by the classical Liouville transformation. We write the Poisson bracket 
and the flows of the CH hierarchy in terms of the potential of the string spectral 
problem. Section 4 presents standard facts about the spectral theory of the string. 
Jacobi ellipsoidal coordinates are discussed in Section 5. The Poisson bracket for 
the Weyl function is computed in Section 6. Action-angle coordinates are intro- 
duced in Section 7. There also CH flows are linearized. Section 8 contains examples 
of the computation of the Weyl function for one and two peakon solutions. Section 
9 introduces generalized action-angle coordinates which include the classical ones 
as a particular case. Section 10 contains another family of canonical coordinates. 
Both families cover all possible spectra which arise from various mixed boundary 
conditions. Section 11 presents the derivation of the trace formulas. 

The author would like to thank Henry McKean for many constructive remarks. 
Many thanks to Mike Frazier for reading the manuscript and suggesting a lot of 
improvements. 



2. The spectral problem 



We consider the Camassa-Holm equation 





— oo 



where t > and — oo < x < oo. The integral operator R defined by 




— oo 
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is an inverse of L = 1 — D 2 . Introducing m — L[v] we specify the phase space M. 
as all smooth positive functions with such decay at infinity that 



oo 

J m(x)e^ dx < oo. 



We also need Mo G M. the subspace of functions which vanish far enough to the 
left. M. and A4o are invariant under the CH flow, [9]. 

The first four Hamiltonians, integrals of motion, are listed in 1.2. The first 

+oo 

Hamiltonian Hi = f v <ir=momentum produces the flow of translation: 

— oo 

= m* + {m,Hi} = rrf + Dm, (2.1) 

where the bracket is 

{A, B\ = [ ^(mD + Dm)^-dx. (2.2) 
J dm dm 

— OO 

The CH Hamiltonian H 2 = \ f v 2 + Dv 2 dx =energy is the second in the infinite 
series of commuting integrals of motion. The equation can be written as 

d 

= m* + {m, H 2 } = m* + (mD + Dm) v = (— + 2(Dv) + vD)(l - D 2 )v. (2.3) 

C6 

All other integrals produce higher flows of the hierarchy. 
The CH equation is a compatibility condition between 

D 2 f - \ f + Xmf = (2.4) 

and 

r = -(v + ±) Df+ l -{Dv)h 



2XJ J 2 

i.e., (D 2 f)* = D 2 (f) is the same as 2.3. The isospectrality of 2.4 is a key to the 
integrability of the CH dynamics. 

3. The Liouville correspondence 
The standard Liouville's transformation 

x - iU) = 2 tanhx/2, f(x) ^ /(0 - J " " ' 



coshx/2 



converts 2.4 into the string spectral problem 



f" + \gf = 0, (3.1) 
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with g(£) = m(x) cosh 4 x/2 and — 2 < £ < +2. The transformation changes the 
length element by the rule 

e 2 

dx — ► d£ = dxj(£), with j-(f) = i _ L.. 

Under the assumptions made, the phase space of the string is regular, i.e., its mass 
is finite: 

+2 +oo 



J 9{£)d£= J tti(x) cosh 2 x/2 dx < oo. 



The transformation reduces the problem 2.4 with two singular ends to the reg- 
ular string problem on the finite interval. It is necessary for introduction of the 
ellipsoidal coordinates but first we consider what is happening on the phase space. 
The functions m(x) and v(x) are related by m = L[v] and v = R[m\. Introducing 
= v(x) cosh 4 x/2, we have the relations 

g(0 = C[k}(0 = - 4 ^' 2 - 2 J J") - WJJ' ~ ^"J\ 

and 

+2 



-2 



with 



i-v/2 \ 2 _j_ forr7>f 



The functions k(£) and g(£) are new coordinates on the phase space. Best of all, 
the relation between the functions m, v and g, k can be expressed by the diagram: 

m(x) : m = L[v] < ► v(x) : v = R[m] 



g(£) = m(x) cosh 4 | J | fc(f ) = v(x) cosh 4 



9(0 ■■ 9 = C[k] 



<— fc(£) : fc = ^b] 

15 



Now we want to express the first two flows 2.1 and 2.3 of the CH hierarchy in 
:ms of the new coordinate, the density function g. The conserved quantities are 



terms 

+2 

Wi = J gjd£, 

-2 

+2 



H 2 = X - J g n[g]J 3 dZ, etc. 

-2 

The bracket 2.2 takes the form 6 

+2 

{A, B} = J*£ (Jgdj- 1 + J-'dgj) ^ d£. (3.2) 

-2 

Computing the variation 5Hi = J ^ 5g J d£ we write the translation flow 2.1 as 

= g - + {g,H 1 } = g' + J- 1 d{gJ 2 ). 
The self-adjointness of 1Z acting on the functions / and g of the ^-variable, 

+2 +2 

J K[g] fJ z di = J gn[f]J 3 d£ 

-2 -2 

implies the formula for the variation of the basic Hamiltonian 

SH 2 = \ J 5gn[g] J*di + \ J 9 n[5g] = J 5g U[g] J^. 
Therefore, for the CH flow 2.3, we have 

= g* + {g, H 2 } = g' + (gjdj' 1 + J- l dgJ)kj\ 
Evidently the Liouville's transformation complicates the formulas for the flow. 

4. The spectral theory of the string 

In this section we introduce the so-called Weyl function, the main tool in con- 
struction of ellipsoidal coordinates. We explain its relation to the spectral theory 
of the string. The direct and inverse spectral theory of the string with nonnegative 
mass was constructed by M.G. Krein in the 1950's and presented in [15], see also 
[13, 10]. To formulate the results we need some facts of function theory. These 
can be found in [14]. 



6 <9 stands for the derivative in £ variable. 
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The function of a complex variable F(X) belongs to the class (R) and is called 
an i?-function if it is holomorphic on C — K. with F(X) = F(X) and ^sF(X) has the 
same sign as 3 A. Any such function can be represented in the form 

+00 

F(X) = a + [3X+ J _ da(C) (3A + 0), (4.1) 

—00 

where (5 > and a is real constant. The positive measure da is such that 



+00 




< 00. 



—00 



Denote by (R) the class (R) adjoined with the function identically equal to infinity. 
Evidently if F(X) € (R) then -1/F{\) G (R). 

Let (i?i) be a subclass of functions which admit an absolutely convergent integral 
representation 

+00 

F(X)=a+J^- (3A^0), (4.2) 

—00 

where a is a real constant and da is a positive measure. 

We introduce two solutions </?(£, A) and A) of the eigenvalue problem 3.1 
with the boundary conditions 

V?(-2,A) =1 V(-2,A) = 

^'_(-2,A) =0 ^(-2,A) = 1. 

The functions ip(2,X) and ^(2, A) can be written in the form 

^) = n(l-£). «2,A) = 4n(l-£). (4.3) 

The roots interlace each other 

< iii < Ai < /i 2 < A 2 < 

Given some function iV(A) which belongs to the function class (R). This function 
describes the boundary condition at the right end of the string. For example, if 
the right end of the string is fixed then N(X) = and vice-versa. The function 
il N was introduced by H. Weyl, [28], in his study of the Sturm-Liouville problem. 
It is defined such that the corresponding Weyl solution x — —i> + satisfies 
boundary condition 

X ' + (2)iV(A) + x(2)=0 for all A. 
The Weyl function is the main ingredient of the spectral theorem. 
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To formulate the spectral theorem we consider the measure dM(£) = g(£,)d£, 

(2) (2) 

where g{£) > is the function from 3.1 and L M = L M [—2, 2] the set of im- 
measurable complex-valued functions /(£) such that 



Lw , |/(0|W(0<oo. 

Evidently, L M is a Hilbert space. For any nondecreasing function u(A) given on the 
interval (— oo, +oo), we denote by L a = L a '(— oo,+oo) the set of a-measurable 
functions F(X) such that 



\F\\l= J \F(X)\ 2 da(X) < oo. 



The string with a fixed left end is denoted by Sq and the string with a free left 
end is denoted by <Si. Consider the string Si first. The nondecreasing function 
a is called a spectral function of a string Si if the mapping U : / — > F, where 
fe L$ [-2, +2], and 

2 

F(X) = J /(Ov(e,A)dM(0, 

-2 

isometrically carries the space L M [—2, +2] into the Hilbert space L a . Accordingly 
the Parseval equation 

+oo 2 

\F(X)\ 2 da(X) = I \f(0\ 2 dM(0 



holds. The spectral function a of the string is said to be orthogonal if the mapping 
U : / — > F maps onto the entire space L^. 

We say that the string has heavy endpoints if £ = —2/ + 2 are points of in- 
crease/decrease of the function M(£). 

Theorem 4.1. The Spectral Theorem. Suppose the string Si has heavy 
endpoints and does not carry a concentrated mass at the right endpoint. Then the 
formula 
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defines an Ri-function for any choice of N(X) G (R). Any such function admits 
the integral representation \.2 with a = 0: 



tt N (X) 



d(x N (C) 



These formulas establish a one-to-one correspondence N <-> a between the class 
(R) of functions N(X) and the set of all spectral functions er(A) of the string Si. 
The spectral function crjv(A) will be orthogonal if and only if the function N(X) 
corresponding to it degenerates into a real constant, possibly infinity. 

To state some facts about strings we need two additional subclasses of (R). 

We say that the function F(X) £ (R) belongs to the class (S) if it is holomorphic 
on (—00, 0] and > there. All such functions have an absolutely convergent 
integral representation 

' + °° da(C) 



F W = <* + J Q ( 4 - 4 ) 

where a > and da is a positive measure. Evidently (S) C (R\). Denote by (S) 
the class (S) adjoined with the function identically equal to infinity. 

Similar, the function F(X) £ (R) belongs to the class if it is holomorphic 

on (—00, 0] and < there. All such functions have an integral representation 

F(X) = a + (3X + jT 00 - dp(C), (4.5) 

where a < 0, j3 > and dp is a positive and such that 

Denote by (5* -1 ) the class (S^ 1 ) adjoined with the function identically equal to 
infinity. It can be proved that F(X) E (S) if and only if -1/F(A) G (S" 1 ). 

Remark 4.2. Let the right end of the string is fixed (iV(A) = 0). For such a string 
Si the Weyl function fi £ (S) and according to the general formula 4.4 



m _ ^(2,A) _ v a k 
fi ° (A) -^A)-^^A- 



k= 

Remark 4.3. If the left end of the string Si contains an interval free of mass then 
the Weyl function takes the form 4.4 with a = I where the constant / is the length 
of the interval. Therefore, for the string Si with fixed right end, we have 

00 

^o(A)=/ + ^^-. (4.6) 
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From formula 4.3 we have 



no( 0) = 4?4 = 4. 
1 ' <p(2,0) 



Combining this with formula 4.6 we obtain 



oo 



l + m =, ,4,) 



Remark 4.4. Given /z's which are spectra of the string S±, what I and what sequence 
of cr's can appear in formula 4.6? For any sequence of a such that 



oo 



— < oo, 
Atfc 



fe=i 



one has a function f2 £ given by formula 4.6. If / is such that fio(O) = 4 (what 
is equivalent to 4.7), then to such Qq corresponds a unique regular string Si with 
heavy right end and the interval at the left of length / free of mass. 

Remark 4.5. For the string So with fixed left end, the set of all Weyl functions is 
given by the formula 

E N = - — . 

If the right end is also fixed then E (X) E (S^ 1 ) and according to general theory 

*w-«sa35-^-s + g(3cb-5;)* (4 ' 8) 

The poles of E (\) are zeros of ip(2, A) and as we will see they do not move under 
the CH flow. The sum 



oo 
k=l 



Pk ^ 

< OO 



if and only if / > 0. For such a string <S , the formula 4.8 becomes 

1 +oo 

k=l K 

Evidently 



I ^ A fe 4 

fc=i K 
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Remark 4.6. Given the sequence of A's which is the spectrum of the string <S , 
what / and sequence of p can appear in formula 4.9? For any sequence of p such 
that 

oo 

one has a function E defined by formula 4.9. If I is such that E (0) = — | (what 
is equivalent to 4.10), then E (— x) < E (0) < for any x > 0. Therefore, 
Eq G (S^ 1 ). The function 

belongs to S. The function VLq has the form 4.6 and satisfies the condition fio(0) = 
4. According to Remark 4 there exists a unique corresponding regular string Si 
with heavy right end and the interval at the left free of mass. 

5. Jacobi ellipsoidal coordinates 

This section contains information about Jacobi coordinates in finite and infinite 
dimensions, see [12] lecture 26. The results are used to describe the range of the 
canonical variables constructed in the subsequent sections. The present discussion 
(and notations) are independent from all other sections. 

First we explain the construction of Jacobi coordinates in finite dimensions. 
Consider the function 

N 

T(A) = £-% u k >0. 

Assume that the poles 7i < 72 < • • • < 7jv are ordered and fixed. The residues 
1/1, . . . , i/jv are considered as variables and they fill up R+. Pick C ^ and consider 
points where T(A) = C . We have N real roots Xk = Xk(C), k — 1, . . . , N. These 
are Jacobi coordinates. If C > then, 

Xi < 7i < X2 < ■ ■ ■ < Xn < In] (5.1) 

and if C < 0, then 

7i < Xi < 72 < • • • < In < Xn- (5.2) 
The converse is also true. Pick some C and consider the points Xk, k — 1, . . . , N 
which satisfy inequalities 5.1 or 5.2 correspondingly. There exist a unique sequence 
of ui, . . . , i/jv from corresponding to the sequence of x's- 
The following formula 7 is due to Boole, [6]: 

N 

Cj2"?k-Xk = s , C^O, (5.3) 

fc=i 



7 A. Volberg pointed this to me. 
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where 

s n = Y^ Vkl ki n = 0, 1, . . . . 

k 

In fact Boole's formula is the first in the infinite sequence of identities. The second 
that we found is 

N 

k=i 

One can continue and obtain the formulas for 

N 

k=i 

in terms of the moments sq, . . . , s p _i for any integer p. We present two proofs of 
Boole's formula. The first is classical and can be found in [19]. 
The identity 

C-T(A) = C ^' Xfc j 

n( A - ik) 

implies 

c i[(x - 7*) - e * ri( A - 7 P ) = c n( A - x*)- 

fe fc PT^fe k 

The coefficients of A^ of both sides match. Matching coefficients of A^ -1 we obtain 
5.3. 

The second proof is a little more complicated, but it can be easily generalized. 
Note that near infinity, 

T(A) -"T p"a^ •••• 



Consider the contour integral 



lm J 



2ni J T(A) - C 

over a circle of sufficiently large radius. The poles and zeros of the function T — C 
produce a contribution to the integral equal to 

^2ik - Xk- 

From another side, near infinity 

^ ( A ) \-2„ in \-3/o„ in n 2 /n2\ 



T(A) - C 
Thus the integral is equal to 

so/C 

which is Boole's formula. 
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-\- 2 s /C - \- 6 {2 Sl /C - s 2 /C 2 ) - 



The second and all higher identities can be proved by contour integration with 
A replaced by A p , p > 1. 

Consider now the infinite-dimensional case. Let 

T(A) ^ 



Ik — A ' 

k=l m 

where 7^ — > +00 as k — > 00 and s = ^z/ fc < 00. For any C we have a sequence 
of x's which interlace the sequence of A's similarly to 5.1 or 5.2. Pick some C^O. 
Can any sequence of interlacing A's arise from some sequence of z/'s similarly 
to the finite-dimensional case? In fact, there is an asymptotic condition on \' s 
which follows from the next result. 

Theorem 5.1. [17] The x ?s are coordinates, i.e. they arise from some C ^ 
and some summable sequence of v s if and only if 

+00 

It* - x*i < °°- 

k=i 

If this holds then 

+00 

C^2(lk- Xk) = so- 
k=i 

This result can be obtained from the rational case in the limit when the number 
of poles tends to infinity. 

Unfortunately the condition so < 00 does not hold in most cases we consider. 
For the string with light left end (i.e. without concentrated mass) we have only 

+00 

— < 00. 

k=i^ 

Theorem 5.2. [23] Pick some C > 0. The x' s « r e ellipsoidal coordinates, i.e. 
they arise from some sequence of u s with s_i < 00 if and only if 



+00 1 



< 00. 



k=i 

If C < , then the above condition should be replaced by 

Proof. Without loss of generality we assume that C = T(0) > 0. Then, 

T(A)-T(0) = expF(A) 
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1 for 7 fc < A < Xk+i, k = 0, 1, . . . ; 

for Xk < A < 7 fc , fc = 1, 2, ... . 



where F e Using 4.1 we have 

T(A) - T(0) = Aexp (| MC) ^!^ ) 

where 

/o(A) = 

We put formally 70 = —00. Writing f = 1 — /, we have 

Now pass to the limit with A = —x — > —00 along the real axis. Then the limit 
exists. This implies 

00 

J^log — < 00, 

k=2 Xk 

which is equivalent to the condition of the theorem. Moreover, 
T(A) - T(0) = -T(0) exp (- £°° ^f^j . 

Conversely, pick C > and the sequence of x's which interlace the sequence of 

7's: 

= Xi < 7i < X2 < ■ ■ ■ , 
and satisfy the condition of the theorem. Define the function T(A) by the identity 

where 

00 

/(0 = E Wi(c). 

k=l 

Evidently, T(A) maps the upper half-plane into itself. It also vanishes at infinity 
and has poles at 7^. Thus, applying 4.1 again 

T(A) = a + /3A + V 1+7fc x A i Uk 2 . 

lk-X 1 + 
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Now pass to the limit with A — ► — oo along the real axis. The existence of the 
limit implies that f3 — and 



oo 



t.1 1 + ->i 

Therefore, since T(A) vanishes at infinity, 

oo 

/I Ik \ 



< + OO. 



k 

00 

i a — 

k 

00 

Vk 



E\ 00 
ikVk \ Vk 
I + 7, 2 Z - 



When C > we have Xk < Ik for all k, and the condition 

+00 , s. 

is equivalent 

+00 

nlk ^ 
— < 00. 

k=i Xk 

This product can be expressed in the spirit of Bool's formula as 

~ Ik C 



k=l * k C ~ 

To prove the formula we consider an algebraic situation with A^ < 00. We 
assume C > T(0) = s_i and this implies that xi > 0- The function 

log(T(A) - C) 

is single valued on the plane cut along the segments [xk,lk], k — 1, . . . , N. Con- 
sider the contour integral 



J-Jl xiogm X)-C) d A 



over a circle of sufficiently large radius. Evaluating this integral when the radius 
tends to infinity, we have 

log(-C). 
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Evaluating it in the finite part of the plane we have to account a contribution from 
the origin and small circles surrounding the cuts. At the origin 

^ J I log(T(A) -C)d\ = log^! - C), 

and over a circle surrounding [xfc,7fc] 

Thus 

JV 

log(-C) = log(s_i - C) + (log7 fe - logxfe) • 

fc=i 

The result follows by exponentiating this equation. The formula extends to all 
other values of C and to functions T with an arbitrary number of poles. We do 
not dwell on this. 

6. POISSON BRACKET FOR THE WEYL FUNCTION 

The goal of the present section is to compute the Poisson bracket 3.2 for the 
Weyl function. 

Theorem 6.1. Let N(() = N be a real constant, possibly infinity. Then, 

{n No (\), n No (n)} = (fiAr (A) - n No (v)f ■ (6.1) 

This fact is central for all further discussion. The proof follows the steps devel- 
oped in [26] for the Dirac operator. We start with two auxiliary lemmas. 



Lemma 6.2. The gradient o/Qjv (A) is 

dn No (\) 

where \ — ~~ 4> + Qn <P- 



Proof. The proof is now standard, see [26] . 

Let / = /(£, A) and s = s(£, /i) be arbitrary solutions of the string spectral 
problem 

f" + \gf = 0, s" + figs = 0. 

Lemma 6.3. Let J = J~(£) be an arbitrary function. Then 

s 2 9 Jd £ - fgjd £ = -^—d [f's - fsf . 

J J jl — A 
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This formula can be verified directly. 

Proof of Theorem 6.1. Using formula 3.2 for the Poisson bracket, Lemma 6.2 
and 6.3 and integration by parts, we compute 



HO 

V J X 2 (A) (JgdJ- 1 + J~ l dgj) X 2 (/i) di 

-\n f\\^)gJd^-x 2 {\)gJd^di 



— WWxM - x(A)xV)] 2 1+2 



2' 



x'Wx(») - x(A)xV) 



x(A)x(/^)- 



// — A 

To compute the value at the point £ = 2 we write 

Y(A) xV) 

.X(A) 

Note that x'-^o + X = & t£ = 2. Therefore, at this point 

XXA) _ xV) _ » r 
X(A) x(A*) °' 
and the contribution from the upper limit vanishes. 

At the lower limit £ = — 2, using x — + ^n ^P we have 

x(x) = n No (x), x'(A) = -i. 

Finally, 

A/ ' W(\)xM - x(A)xV)] 2 U = t~~~ (SWA) - n N Mf 



ii — X X — ii 

The proof is finished. 

If one changes the spectral parameter by the rule 

\' = -\, 1*' = —, 

then the formula of the theorem becomes 

{n No (x'),n No (,')} = ^y)-W))\ (62) 

A — Li 

This is the standard Atiyah-Hitchin bracket, [3], in the form found by Faybusovich 
and Gehtman, [11]. Note that the function fijv (A') maps the upper half-plane into 
itself. Due to the invariance of the AH bracket under linear-fractional transfor- 
mations, [26], the bracket for the function E No (X') is given by the same formula 
6.2. 
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7. Action- angle variables. Dynamics of the Weyl function 
The main result of this section is 
Theorem 7.1. The variables 

I n = —r-, e k = logy?(2, A fc ) 

are "action-angle" variables 

{I k ,I n } =0, 

{#fc, #n} = 0, 
{9k, In} — 



We start the proof with two auxiliary lemmas. In section 4 we introduced the 
representation 4.8 for the function E (X): 



Pk- 



k=i 

Introducing 



K — Pk — T2, (7-1) 



we have the identity 



Afc — A X' k — A' Afc 
Therefore for the function E (X') we obtain 



Pk _ P' k + Pk^ 



i 00 / 

^(A') = -t + Ev^T7- ( 7 - 3 ) 



fc=i K 

First, we compute the bracket between the the parameters X' k and p' k entering 
into formula 7.3. 

Lemma 7.2. The bracket 6.2 in X' and p' coordinates has the form 

{Pk',p'J = 3r4r(i-tf)> (7-4) 

A n ~ A k 

{p' k ,X' n } = p> k 6 n k , (7.5) 

{X' k ,X' n } = 0. (7.6) 

The proof is identical to the proof of Theorem 2 in [25] . 

Now we can compute the bracket between the parameters Afc and p k entering 
into formula 4.8. 
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Lemma 7.3. The bracket 6.1 in X and p coordinates has the form 

f 1 ^X k X n p k p n /„ „\ 

iPk,pn} = — : r , (7.7) 

{p k ,X n } = Xt Pk 5 n k , (7.8) 

{X k ,X n } = 0. (7.9) 

The proof follows from the results of the previous lemma. 
Proof of Theorem 7.1. The first identity follows trivially from 7.9. Commuta- 
tivity of angles is more complicated. Using 4.3 



R,U = ^{iogfi-M,io g fi-M} 



_ sr^ 1 | A n X m ^ 

~ ^1 _ ^) (l - ^ ) A*fc'^fc' 

= 52 7 ; \ / " "1 \ {{Hk, Pk'} + {Pk', Pk}) 

= 0. 

The last sum vanishes due to skew symmetry of the bracket. 

To prove the last relation we use the identity <p(2, X k ) = ^'(2, X k )p k . We have 

{9 k ,I n } = {log^'(Afc) +logp fe ,-^} 

- 1 -{Pk,Xn} = ^Xl Pk 5 n k =6 n k . m 



PkX 2 n p k xi 

H. McKean, [20] expressed the conserved quantities Tii and H2 in terms of the 
Dirichlet spectrum 

7° 7 +°° 

Hi = mdx= g Jd^ = s ^X^ 1 , 
—00 — 2 

+°° + 2 +00 

H 2 = l -J mR[m}dx= l - J gTZ[g]J z d£ = ^ A; 2 . 



These are obtained by matching coefficients of two expansions near A = for 
the solution ■0. One expansion is obtained from formula 4.3; the another form a 
Neumann series 8 . 



These and more general trace formulas will be obtained in Section 11. 
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The evolution of residues under the translation flow 

pl = {pk^Hi} = ^{p/cA; 1 } = -X k 2 {p k ,\ k } = -p k , 

n 

i.e. 

Pk(t) = p fc (0)e~*. 
Therefore, under the translation flow 



E (\) = - 




The evolution of residues under the CH flow obeys 

1 x 1 

Pl = {Pk, ^2} = 4 Yl = ~ 2 Afc ^ = ~PkPk, 

n 

where p k = Therefore, 

Pk(t) = Pk(0)e-^\ 

and 

8. Examples. 1 and 2 peakon solutions 

A remarkable class of solutions of the CH equation are peakon-antipeakon so- 
lutions of the form v(x,t) = J2 n p n e~^ x ~ qn ^ The parameters q n and p n satisfy the 
Hamiltonian flow 

Qn = dH/dp n , p' n = -dH/dq n 

with Hamiltonian H — | YliPiPj e ~^ qx ~ q ^ an d the classical Poisson bracket. 

It is instructive to compute the evolution of Weyl functions for the simplest 
one and two peakons solutions. The equations of motion of two peakons were 
already integrated in [7], but the evolution of the Weyl function computed here is 
an illustration of results of the previous section. Since everything is algebraic we 
allow p n to be of both signs. 

Example 8.1. One peakon. 

The solution has the form v(x, t) = p(t)e~^ x ~ q( - t ^, with pit) = po and q(t) = qo+pt. 
The peakon (p > 0) travels to the right while the anti-peakon (p < 0) travels to 
the left. 

The solutions ip and ip take the form 

(p{2, A) = -AmiZi + l, 
■0(2, A) = -AmM + 4. 
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Using miWi = 8p for the Weyl function of the string So with fixed right end 
(N = 0) we obtain 

F (\\ V(2 ' X) - 1 4- Pit) 



■0(2, A) l i-A' 



2p 

with 

i i P -?o 

p(t) = — - = _ e - q = - — e~ pt . (8.1) 
yK J mi/g 8p 8p K ' 

Example 8.2. Two peakons. 

The solution is of the form 

v(x,t) = Pl {t)e-\ x - q ^\ + p 2 (t)e-\ x - q ^\ . 
The two-particle Hamiltonian is 

where p n is the asymptotic velocity of the n-th particle. The equations of motion 

OH OH 
ql = —= Pl +p 2 e-\^\, pl = - —= p lP2 sign( qi -q 2 )e-^-^, 
Opi Oqi 

OH OH 
q' = — = p 2 + p!e~\ qi ~ q2 \ p* 2 = -— = -pip 2 sign(gi - q 2 )e~ lqi ~ q ' 21 , 
0p 2 0q 2 

are integrated by introducing the new variables 

P = pi +p 2j Q = g 1 +g 2) 

P =Pi-P2, q = qi~q2- 

This reduces reduces the equations of motion to 

P* = 0, Q 9 = P(l + e-\ q \), 

p" = ~{P 2 -p 2 ) signge" 1 " 1 , q* = p (1 - e~M). 

Now the Hamiltonian is 

H = ±[P 2 +p 2 + (P 2 -p 2 )e-\ q \}. 

Consider the case 9 p 1 > p 2 > 0. This is the pure peakon case. The law of 
conservation of energy takes the form 

AH -P 2 - p 2 _ k | 
P 2 -p 2 ~ 6 ' 



9 The case p 1 = p 2 is excluded by the equations of motion. 
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with left hand side less then 1 for all times. This implies that q = qi — g 2 < for 
all times. The solution is 

„ 1 - P eM 
p(t) = A 



q(t) = log 



A 2 (3e At 



(Pi + P2p eAt ) iP2 + Pi(3e M ) ' 



r>u\ nm\ , d + i Pi+P2P eM . , P1+PP2 

Q(t) = e(0) + ^- log STfA^ + log ST^ i 

with 

The solutions ip and ip are 

ip(2,X) = A 2 /i/ 2 ^i^2 — A(m 1 / 1 + miZ 2 + 7?7. 2 / 2 ) + 1, (8.2) 
ip(2,\) = \ 2 m 1 m 2 lohk - A(mi/ / 2 + miWi + m 2 lol2 + m^hh) + 4. (8.3) 

Note that 

777.1777.2/0^2 = 8P 2 — 16i/ = I6pip 2 

and 

777i/ / 2 + 777i/ /l + ?™ 2 / / 2 + = 8P = 8(p x + p 2 ). 

Therefore, for the poles (A : if)(\) = 0) we obtain 

Ai = — A 2 = — . 

2pi 2p 2 



Thus: 



with 



7?(2,A) 1 Pl (t) p 2 (t) 

tio{A) = -77TTT = ~T + "1 r + 



and 



^(2, A) " k 4-A 4- A 



* _ /J P 1 + ^ 2 ' <8 ' 6) 
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Peakon-antipeakon. Now to the case p 1 > > p 2 . The integration of the 
equations of motion is identical to the pure peakon case. The solution is 

1 + (e At 



p(t) = A 



q(t) = log 



A 2 (e At 



{Pi(e At ~P 2 )(Pi -V2& M Y 



where 



Thus, 



with 



and 



Q(t) = Q(0) + Pt-log Pl _ p2 _ C f +log P2 _ P ^f : 

Pi - P2C P2 - PiC 

p(0) + A 
ip(2,X) 1 Pl (t) p 2 {t) 

= p 2 (t) = 



e -O(0)p iC _p 2 



C Pi - P2C 

The formulas for residues are obtained in the same way as in the pure peakon case. 
At the moment of the collision the mass of the first particle tends to infinity and 

si = pi(t)+ /02(f) = 0- 

9. Generalized action-angle variables 

The action-angle variables introduced in section 7 are a particular case of the 
following construction. According to formula 4.6, the function Jlo(A) has the form 

All possible sequences of /i's are described in Theorem 1.2. 

Pick some constant C which does not exceed /, and find the points where fi (A) = 
C (see Figure 1). We have the sequence of points = Afc(C), k — 1, 2, . . . which 
depend on the constant C. They interlace the poles of ^o(A): 

< Hi < Ai < fi 2 < A 2 < • • • • (9.1) 

As it was explained in section 5, these are Jacobi ellipsoidal coordinates. Section 
7 describes the case when C — 0. Given /i's, what is the range of all possible 
sequences of A(C)'s? 
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The Nazarov-Yuditskii theorem implies 



Conversely, let the sequence of A's satisfy 9.1 and the above condition. Then again 
by the Nazarov-Yuditskii theorem there exists a function 



such that T(Afc) = —1, k — 1,2, Note that T(0) > depends only on the 

given sequence of //'s and A's. Pick I, < I < 4 and define 

w )=l+ (±^m. ( , 2) 

By construction 

^o(A fe ) = /-^^ = C, fc = l,2,... 

the range of C being the interval — 4/T(0) < C < 4. The function fi(A) has the 
form 4.6 and satisfies condition 4.7. Therefore, according to remark 4.4, for any 
choice of /, or equivalently C, there exists a string Si with the function fi given 
by formula 9.2. 

Remark 9.1. For C = 0, there exists a string 5i with / > and corresponding 
Dirichlet spectrum coinciding with the A's. 

Now we want to construct variables canonically conjugate to A(C). 

Theorem 9.2. The variables 

In{c) = ~\Jcy 6k{c) = logv?(2, Xk{c)) 

are canonically conjugate variables 

{I k (C),I n (C)} =0, 

{0 k (c),e n (c)} =o, 

{9 k (C),I n (C)} 



To prove the theorem we introduce the function Eq(X): 

1 



n -c 
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It has the representation 

<p(X) <p(X) 



E C (X) 



V(A) - C<p(X) V C (A) 

The points Xk(C) are the roots of the equation ip c (X) = 0. The function ip c (X) can 
be written as 

nx) ■ - 1T ' A 



(4-on 1- 

k=i ^ 



A fc (C) 



Evidently 4 — C > 0, since C < I < 4. The function -E'g(A) maps the upper 
half-plane into itself and can be written as: 

BS(A) ._^ + g_|M_. (, 3 ) 

One can change the variable A — > A' = — ^ and get the formula 

After that the proof of the theorem follows familiar steps. The Poisson bracket 
for the function Eq(X') is given by formula 6.2. For the variables p'(C) and X'(C) 
from 9.4, an analog of Lemma 7.2 holds. Then similarly to Lemma 7.3 we can 
compute the bracket between the variables p(C) and A(C) entering into formula 
9.3. Finally, the proof of the canonical relations is identical to the proof of Theorem 
7.1. ■ 

We will construct a family of Hamiltonian flows which will be linearized in the 
variables 1(C) and 9(C) in section 11. 

10. Second family of canonical coordinates 

The arguments of the previous section prompt the construction of another family 
of canonical coordinates J and r. As it was explained in the introduction these 
variables are complimentary to / and 9. We assume that the initial data belong 
to the smaller phase space M.q. 

According to formula 4.8, the function Eq(X) has the form 

<p(2,\) 1^ p„ 



ip(2, A) I ^A fc -A 

The range of all possible sequences of A's is described in Theorem 1.1. 

Pick some F such that — j < F and consider the points where -Eo(A) = F (see 
Figure 2). We have the sequence of points Hk = l^k(F), k = 1,2, ... . interlacing 
the poles of E (X): 

fi! < Ai < ^2 < A 2 < . . . . (10.1) 
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Given A's, how can the range of all possible sequences of //(.F)'s be described? 
The Nazarov-Yuditskii theorem for this case implies 



k=i 



Conversely, suppose the sequence of /i's satisfies 10.1 and the above condition. 
Then again, by the Nazarov-Yuditskii theorem, there exists a function 

oo o oo o 

such that T(/t fe ) = 1, k = 1, 2, Note that Ai > and T(0) > depends only 

on the chosen sequence of /i's and A's. Pick I, < / < 4 and define 

F(\)- 1 i ( 4 ~0 T ( A ) ri0 2 ^ 

Eo(A) -"T + 4/T(0) • (1 °- 2) 

By construction 

em= -j + M) =f ' k=i ' 2 -- 

For the range of F we have three options. 

• If /ii > then T(0) < 1 and F satisfies —1/4 < F < +oo. 

• If iii < then T(0) > 1 and F satisfies -oo < F < -1/4. 

• If /i! = then T(0) = 1 and F = -1/4. 

The function E (\) has the form 4.9 and satisfies condition 4.10. Therefore, ac- 
cording to remark 4.6 for any choice of /, or equivalently F, there exists a string 
So with the function E given by formula 10.2. 

Remark 10.1. Let /xi > and F — 0. Then there exists a string <S with / > and 
corresponding spectrum /i(0)'s. 

The next theorem is an analog of Theorem 9.2. We pick F such that H\(F) ^ 0. 

Theorem 10.2. The variables 

J n (F) = }— , r k (F) = log ^(2, n k (F)) 

are canonically conjugate variables 

{J k (F),J n (F)} = 0, 
{r fe (F),r n (F)} = 0, 
{r fe (F),J n (F)} = 
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To prove the theorem we introduce the function fi^(A): 

1 



(A) 



oo , 

^(A) = (l + 4F)n 1 

k=l ^ 



E (\) - F 
It can be written as 

The points Hk(F) are the roots of the equation y? F (A) = 0. The function (p F can 
be written as: 

A 

M>). 

The function £Iq(\) maps the upper half-plane into itself and has the expansion 

One can change the variable A — > A' = — ^ and get the formula 

After that the proof of the theorem follows familiar steps. The Poisson bracket 
for the function ^Iq(X') is given by formula 6.2. For the variables cr'(F) and fJ,'(F) 
from 10.4 an analog of Lemma 7.2 holds. Then similarly to Lemma 7.3 we can 
compute the bracket between the variables cr(F) and fJ>(F) entering into formula 
10.3. Finally, the proof of the canonical relations is identical to the proof of 
Theorem 7.1. ■ 

We will construct a family of Hamiltonian flows which will be linearized in the 
variables J(F) and t(F) in the next section. 

11. Trace formulas 
Let Afc(a, b) be a spectrum of the boundary value problem 

aip(2,\) + 6^(2, A) = 0. 
In this section we will express the sums 

+oo ^ +oo ^ 

^a^M)' £a|M)' etc; 

in terms of the potential. In order to simplify calculations we work with the 
Sturm-Liouville problem 2.4. We assume that the potential m(x) has compact 
support. The solutions ip(x, A) and i^(x, A) of the Sturm-Liouville problem are 
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the image of the corresponding solutions of the string problem under the Liouville 
transformation. Note 

<p(x,0) = (fo(x) = coshx/2, tp(x,0) = ip (x) = 2e x/2 . 

We also introduce the solution 9(x, A) which corresponds to the solution A) 
specified by the boundary condition at the right end 

0(+2,A)=O, 9' + (+2,X) = l. 

Note that 

0(£,O)=£-2, 9(x,0) = e (x)=2e- x ' 2 . 

At the points of the spectrum Xk(a, 6), the solution A(x, A) = aip(x, A) + bip(x, A) 
tends to zero when x tends +00. At these points, the Wronskian 

W(X) = 6(x, X)DA(x, A) - A(x, X)D9(x, A), 

vanishes. 

The expansion can be obtained near A = 

W(X) = I (a, b) + h{a, b)X + I 2 (a, b)X 2 + ... 
in two different ways. The first arises from the Hadamard formula 

W(X) = (a + 46) f[ (l - A\ \ k = X k (a, b). 

k=i V Xk <> 

The coefficients have the form 

I Q (a,b) = a + 46, 

h(a,b) = -(a + 46) J] A 



/ 2 (a,6) = (a + 46) £ -A 



Another expansion is obtained from the Neumann series for the solutions 

V?(x,A) = ip (x) + (pi(x)X + (p 2 (x)X 2 + . . . ; 
tp(x,\) = i>o(x) + *pi(x)X + ip 2 (x)X 2 + . . . ; 
9(x,X) = 9 (x) + 9 1 (x)X + 9 2 (x)X 2 + .... 

Therefore, 

J (a,6) = 9 DA -D9 A , (11.1) 

7i(a,6) = 9 DA 1 + 9 1 DA -D9 A 1 -D9 1 A , (11.2) 
J 2 (a,6) = 9 DA 2 + 9 1 DA 1 + 9 2 DA -D9 A 2 - 09^-09^0; (11.3) 
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where A k (x) = a<pk(x) + bipk{x). Identifying the corresponding coefficients we 
obtain the first two trace formulas 



+00 



^A t (a,6) 

+ OO ^ 



hja, b) 
' a + 46 ' 

h{a,b) 
a + 46 



2/ 2 (a, 6) 
a + 46 ' 



The coefficients y?fc(x) and ipk{x) are given by the formulas 



cosh — , 
2' 



-2 / d{, sinh— ^— cosh-m(£), 



—00 



(11.4) 
(11.5) 



// X — £ £ — 77 77 

dl; / c/?7 sinh — - — sinh — - — cosh - m(£,)m(r)), 



—00 — oc 



4 

2e x '\ 

X 

-4 I di S inh^^e«/ 2 m(0, 



Also, 



^2(2?) 



Oi 



di] sinh sinh e^ 2 m(^)m(ry). 

2 2 



-00 —00 



2e-*/ 2 , 



+00 

= 4 J ^ s inh^^e-«/ 2 m(0, 



+00 +00 



dt] sinh — - — sinh — - — e~ v ^ 2 m(^)m(r)). 



It is easy to check that formula 11.1 produces Io(a,b) = a + 46. Note that the 
coefficients h(a, 6) are linear in the parameters a and 6: ifc(a, 6) = a/fc(l,0) + 
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Wfc(0, 1). For the first trace formula 11.4 we have 



ah(l,0) 6/ 1 (0,l) 



y 1 _ 

X k (a,b) a + 46 a + 46 

K = l 



From 11.2 with the help of formulas for the coefficients of expansions of the solu- 
tions we have 

f+OO 



/+oo 
dxm(x)(l + e- x ), 
-oo 

/+oo 
dxm(x). 
-oo 



For the second trace formula 11.5 we have 

/ a/ 1 (l,0) + 6/ 1 (0,l)V a/ 2 (l,0) + 6/ 2 (0,l) 



y— 

\2( n 



X 2 k (a, 6) V a + 46 / a + 46 



From 11.3 with the help of the formulas for the coefficients of the expansions we 
have 

/ +oo \ 2 +oo +oo 

/ 2 (1,0) = \ \ J dxm(x) J —2 J dx e~ x m{x)D^ 1 m{x) — J dxm(x)v(x), 



— CO / — oo 

+oo \ 2 +oo 



I 2 (0, 1) = 2 J J dxm(x) J -A J dxm(x)v(x). 

\ — oo / — oo 

The first series of Hamiltonians 7ik{C) corresponds to the case a = —C, b = 
1. Two Hamiltonians 7i\{C) and H.2{C) are given in 1.4. The second series of 
Hamiltonians 1k(F) corresponds to the case a — 1, b = F. The first few also can 
be found in 1.4. 
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